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MULTISCALE FUNCTIONS, SCALE DYNAMICS AND 
APPLICATIONS TO PARTIAL DIFFERENTIAL EQUATIONS 

by 

Jacky Cresson & Frederic Pierret 


Abstract. — Modeling phenomena from experimental data, always begin with a choice of hy¬ 
pothesis on the observed dynamics such as determinism, randomness, derivabihty etc. Depending 
on these choices, different behaviors can be observed. The natural question associated to the 
modeling problem is the following : “With a finite set of data concerning a phenomenon, can 
we recover its underlying nature ? From this problem, we introduce in this paper the definition 
of multi-scale functions, scale calculus and scale dynamics based on the time-scale calculus (see 
m)- These definitions will be illustrated on the multi-scale Okamoto ’s functions. The introduced 
formalism explains why there exists different continuous models associated to an equation with 
different scale regimes whereas the equation is scale invariant. A typical example of such an 
equation, is the Euler-Lagrange equation and particularly the Newton’s equation which will be 
discussed. Notably, we obtain a non-linear diffusion equation via the scale Newton’s equation 
and also the non-linear Schrodinger equation via the scale Newton’s equation. Under special 
assumptions, we recover the classical diffusion equation and the Schrodinger equation. 
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1. Introduction 

This article deals with new mathematical tools to deal with scale phenomena and applica¬ 
tions to partial differential equations. The framework that we have developed can be read 
from a mathematical point of view following each definitions and theorems. However, this 
framework can be seen as a synthesis of different tentative of one of the authors in order 
to deal with scales in geometry and analysis in the context of different physical problems 
(see m 0 El la ini), in particular the scale relativity theory developed by L. Nottale 
[2Sl ED and more generally modelling problems. As a consequence, before coming 
to more mathematical considerations, we picture some important problems in modelling in 
Physics which are underlying our framework. 

Modelling a given phenomenon from experimental data using classical mathematical tools 
always assume, sometimes implicitly, a given framework hypothesis on the real nature of 
the phenomenon which can be also called the texture of reality. As an example, classical 
mechanics is developed using the classical differential calculus to write speed and acceleration 
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of particle and implicitly assuming that the behaviour of these particles can be described 
using smooth curves on a given space. Depending on this framework, different behaviors will 
be predicted or not and will be confronted to reality. However, this assumption about the 
real nature of a phenomenon is in general not so easy to decide and in some sense depends on 
philosophical considerations (positivism, etc) which can not be proved. The classical debate 
between A. Einstein and N. Bohr about the nature of quantum physics is a famous example. 

The previous problem can be handled using a different approach, looking at the way math¬ 
ematical models for a given phenomenon are constructed. Indeed, the framework question is 
in fact related to two different facts which are in general mixed in the literature. In order to 
put in evidence these points, we first remind very roughly the usual way to construct a model 
for a given phenomenon : 

Acquiring experimental data. 

- Computations of relevant quantities (velocity, acceleration, etc). 

- Functional relation between these quantities (at a discrete level). 

- Asymptotic passage to a continuous model under a specific choice of hypothesis. 

- Comparison to reality using numerical simulations 

Putting apart the last step for the moment, we see that the framework assumption has to 
do with the following points : 

- Scale dependence : Experimental data are intrinsically scale dependent via the measure¬ 
ment apparatus which induces scale of observation. This remark is well known but the 
mathematical framework to deal with it is a priori not developed. We return to a more 
pragmatic way to deal with this question by introducing the notion of scale functions 
which are basically inhnite family of discrete functions define on a given time-scale. All 
these objects and notions are based on the time-scale calculus introduced in the late 
88 by S. Hilger (see [3]) in order to unify the classical continuous analysis and discrete 
calculus of finite differences. 

- Asymptotic behaviour : The second remark is that the framework assumption concerns 
the asymptotic behavior (if any) of a scale dependent function. This is clearly something 
which is not discussed in the literature because the asymptotic procedure used is not 
put in evidence and implicitly assumed. The usual way is precisely to begin with an 
asymptotic object which is designed in a given framework and to check and rely on 
observed quantities with a given scale corresponding to the precision of measurement. 
In this paper, we define scale equations and an asymptotic procedure on scale functions 
based on the notion of scale regime which allows us to define a natural asymptotic object 
to a scale dependent one, a scale regime being more or less a stable scale behavior of a 
given scale function. 

As a consequence, the framework assumption or texture of reality is a fluctuating notion. 
It depends on the observed scale regime and the initial scale equation. As an example, 
classical mechanics corresponds to a linear scale regime and a part of quantum mechanics to a 
fractional one. We then recover different continuous objects representing a given phenomenon 
on different scale regime. 

It must be pointed out that this result is not trivial. In general, having a given framework, 
a continuous model is written and must, by construction, cover all the scale of observation. 
Doing so, for a given scale, one must understand how some quantities can be negligible. 

In our setting, a given scale regime will modify the corresponding continuous models, 
making some perturbations terms to appear or disappear. A new understanding of the way 
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some equations can be seen as bifurcation or perturbation of some other is then possible. 

This approach, although natural with respect to the modeling problem, leads to the follow¬ 
ing question : If all the quantities and eqnations are scale dependent does it means 
that no nniversal principle can be derived to described the asymptotic equations 
describing a particular phenomenon ? As Physics is based on the search for universal 
laws or principle, this problem is fundamental. The answer is fortunately no. A universal 
principle is a principle which can organize all scales from a given one. An example of such a 
tentative is given by the scale relativity principle initiated by L. Nottale [2011211122] which 
was, as we already mentioned, the inspiration for part of this work. This principle states 
that the equations of motion correspond to some critical point of a Lagrangian functional at 
all scales or in other words that they keep an Euler-Lagrange form. This sentence can be 
rigorously defined using the notion of scale invariant equations in the formalism of the scale 
calculus that we introduce following our previous works H El El [IDl [12] which is based on 
embedding formalism initiated in mm- We are in particular able to define natural analogue 
of classical notions such as Lagrangian, functional, Hamiltonian, symmetries, first integrals, 
etc. 


Observations 



Discrete models 

Continuous models ^- Scale models 

' Numerical ' 

integrators 



Figure 1. Modeling problem. The red connection means that this is dependent of 
explicit choices. 


In this paper, we define what we call multiscale functions which are the geometrical objects 
underlying all modeling problems and the analysis tools necessary to study the behavior 
of these objects under change of scales. These definitions are connected with our previous 
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attempt to define scale manifold in for the geometric part and the scale calculus introduced 
in mm and further developed in m for the analysis part. 

Our constructions and definitions use the framework of time-scale calculus initiated by 
S. Hilger [T71 [T8] in 1988 and further developed in The time-scale calculus unifies the 
classical and discrete calculus in the sense that one can do, for example, variational, integral 
or differential calculus on a continuous interval of time, a discrete interval of time or a mixing 
between both of them. The formalism developed in this paper, called scale dynamics, is 
exactly the feature which comes to supplement the time-scale calculus. Indeed, the time-scale 
calculus only deal with one scale whereas, our new objects allows dealing with multiple scale. 
Also, contrary the previous work concerning the change of hypothesis, scale dynamics allows 
having multiple hypothesis which are encoded in what we called scale laws. It shows why 
there exist different kind of models of the same equation governing the phenomena observed. 

The paper is organized as follows: 

In Section 2, we define scale functions and multiscale functions using the Okamoto’s func¬ 
tions as an illustration. We also define a topology which gives a way to analyse the structure 
of these functions with the help of symbolic dynamics. In Section 3, we discuss the modeling 
problem. Precisely, we come back to the questions of the introduction in order to show the 
necessity of an analysis of multiscale objects. In Section 4, we introduce the scale calculus 
using the definitions and notations of the time-scale calculus. Notably, we define the notion 
of scale equation, scale laws, scale range, scale regime and scale invariance. In Section 5, we 
introduce the scale dynamics. Using the scale and multiscale version of Okamoto’s functions, 
we show the first implication of scale dynamics on these functions. Then, we provide the 
general transformation formulas to link the scale structure together. Precisely, we give the 
formula to quantify the dynamical effects induce by change of scales over the scale derivatives. 
In Section 6 and Section 7, we define asymptotic differential operators and asymptotic scale 
models which are illustrated with the linear and fractional scale regime. This is done using 
the particular way of scale dynamics to decompose a scale or multiscale function as a “reg¬ 
ular” and “irregular” part. These operators depend on the scale regime and the scale range 
chosen to construct the asymptotic continuous model. In Section 8, we apply the formalism of 
scale dynamics to partial differential equations. Using the scale Newton’s equation we obtain 
the diffusion equation and the Schrodinger equation under a change of variable. We can see 
that the diffusion is governed by the special scale regime chosen which is the fractional scale 
regime. 


2. Multiscale functions 

2.1. Okamoto’s functions. — First we remind the classical “one” scale Okamoto’s 
function describe in |26j . 

Let Fa be defined inductively over [0,1] by iterations fi for i > 0 as follows: fo{x) = x for all 
X G [0,1], every fi is continuous on [0,1], every fi is affine in each subinterval [/c/3*, (/c-|-1)/3*] 
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where € {0,1, 2,... , 3* — 1}, and 
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Given this construction, we denote the limit function Fa defined by Fa 
illustrate the construction with multiple a in Figure [2j 


(la) 

(lb) 
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(l d) 

lim fi. We 
2—>-00 


We have (see |23[ Theorem 1 and Corollary 1, p.2-3]) : 

Theorem 1. — Let a G [0,1] and qq is the unique solution o/54o^ — 27a^ = 1. Okamoto’s 
functions have the following properties : 

1. If a G [oo, 1[, then Fa{x) diverges for almost all x G [0,1], 

2. //o g] 0, l/3[U]l/3, ao[, then Fa{x) = for almost all x G 

3. If a = 1/3, then Fa{x) = 1 for all x G [0,1]. 

This result is based on the work of |26] and m 


The next sections introduce a set of notions including the Okamoto’s functions in a wide 
class of objects called scale functions which encode all the classical construction of fractal 
functions. We also define some useful algebraic manipulations on scale functions called scale 
composition. 


2.2. Scale functions. — The previous construction leads to the following objects : 

- Time-scale : Let a, 6 G M, a < 6. A (discrete and hnite) time-scale on [a, b] denoted by 
T is the data of a finite number of points ti G [a, b]. We denote T = {L}. 

- Discrete function : A discrete function is an element of C'(T,R) where T is a given 
time-scale. 

- PL-Continuous representation of a discrete fnnction : We denote by Fpi the lin¬ 
ear interpolation of the discrete function where PL stands for piecewise linear continuous 
functions. 

In this article, we are concerned with more complex objects depending on scale. In order to 
make precise this dependence, we introduce the notion of scale sequences and scale functions. 

Definition 2 (Scale sequence). — Let a,b G M, a < b. A scale sequence denoted by T 
is the data of a one parameter family of (discrete and finite) time-scale T* G [a, 6] such that 
Ti C Ti+i. 

Definition 3 (Scale function). — LetT be a scale sequence. A scale function is the one 
parameter family of discrete functions Fi G (^(Ti,^) denoted by F such that Fj+i jTi= F). 
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(a) 0 = 4/9 (b) 0 = 1/2 
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(c)a = 2/3 (D)a = 5/6 

Figure 2. Examples of Okamoto’s functions 


Figure 3. Scale sequence for Okamoto’s scale function 


As an example, the construction of any Okamoto’s function induces a scale function. 

A three dimensional representation of scale functions can be obtained as follows : for each 
i € N*, we plot on Figure 0] the graph of the PL-continuous representation of Fi denotes by 
bu 

2.3. Scale elementary action. — In order to define a scale composition, we introduce 
the notion of elementary time-scale which are time-scales reduced to only two points, i.e. 
Teiem IS always of the form Te^em = Any (discrete and finite) time-scale can be 

decomposed in a union of elementary time-scales as follows : 

Let T = {to, • • •, tn}- We denote Teiem,i = {ti, ^i+i} for t = 0,..., n — 1. We have 


n—1 


T = 


U Telem,i 


i=l 


( 2 ) 
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Figure 4. PL-continuous representation of Okamoto’s scale function for a = | 


We denote by Elem the operator acting on (discrete and finite) time-scale producing the 
decomposition in elementary time-scales, i.e. 

Elem(T) = {'^elem,i}i—Q^ 

An elementary discrete function is a discrete function defined on an elementary time-scale. 
These discrete functions are the basic piece to describe a given elementary pattern. Precisely, 
we have : 

Definition 4 (Scale elementary action). — A scale elementary action is the data of an 
operator A : C{Teiem,^) C'(TeZem,A, I^) satisfying the following properties : 

- elem,A is such that Teiem C '^eiem,A cind where the intervals associated with the time-scale 
are such that It , a C It , 

•‘‘-eierrn-^ ^elem 

- For all F G C'(TeZem) 1^)? action of A on F denoted by A @ F with A @ F G 
C{Teiem,A,'^) Satisfies 

i^@F)k^ur.= F- ( 4 ) 

An example of scale elementary action is given by the Okamoto’s construction. 

Definition 5 (Scale elementary Okamoto action). — Let TeZem = {^oTi} o,nd a G 
]0,1[. We denote by Oa and we call the scale elementary Okamoto action the scale elemen¬ 
tary action defined by : 

- The time-scale TeZem,Oa • 

'^elem,Oa = {^oTo,l To,2, }, (5) 

pL 

where to,i = to + i-^, i = 1, 2, /r = - to- 

- The action Oa '■ for F G (^(Teiem)!^) we have Oa{to) = F{to), Oa{ti) = F{ti) and 

Oa{to,i) = F{to) + a(E(ti) - F{to)), Oa{to,2) = F{to) + (1 - a){F{ti) - F{to)). (6) 
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In order to see the concrete action of Oa to produce an elementary pattern, we apply Oa 
on the elementary discrete function iiio : I’d = {0,1} —)■ M defined by Eo{0) = 0 and E{1) = 1. 
We obtain the following picture for Oa @ Eq G C{Ti,M) : 



(A) Eo 


2.4. Scale action on discrete functions. — This Section is devoted to the definition of 
a scale action which is the formal formulation of the classical idea of iterative construction of 
fractal function using a given elementary pattern. The elementary pattern is encoded by the 
scale elementary action. The new figure obtained by an iteration of this pattern on a given 
discrete function will be encoded by the scale action. Precisely, we have : 

Definition 6 (Scale action). — Let A be an elementary scale action. Let E E (^(T,]^) 
where T is an arbitrary (discrete and finite) time-scale. 

- Let Elem(T) = {Tj}j=o,...,n-i- For each Tj E Elem(T), as F It^E C(Tj,M) is an elemen¬ 
tary function, the discrete function A{F ItJ E is well defined. We denote by 

the time-scale 

n—1 

TTa = U (7) 

i=0 

The scale action induced by A on F denoted by A @ F £ C'(TA,ffi) is defined by 

[A @ F] A @ [F \t^] (8) 

for i = 0,... ,n — 1. 

The scale action associated to an elementary scale action A allows us to make composition 
of a given operator on a discrete function. Let A be an elementary action. We denote by A®'^ 
the operator defined over discrete function given by the composition of the action oi An times. 

For example, let us consider the iterative action of the scale Okamoto elementary action 
on the discrete function Eq. We denote by Ei = 0®®(Eo) E C'(Tj,R) where Tj = Tj_i^o with 
To = {0,1}. The discrete function E 2 looks like 
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4/9 


2/9 


5/9 


2/3 


8/9 


7/9 


1/9 


1/3 



1/9 2/9 1/3 4/9 5/9 2/3 7/9 8/9 


Figure 5. E 2 = Oa @ Ei = Of^iEo) 


2.5. Multiscale functions. — In this Section we introduce a toy model to explain the 
difficulties related to the modeling of complex phenomenon from a finite set of data. More 
complex constructions can be made and we discuss some of them in the next Section. 

Definition 7 (Multiscale functions). — Let A = {Ai} be a sequence (potentially infinite) 
of scale elementary actions and N = (A^i,... ,Nn) such that Ni € N*, i = 1,... ,n be the 
complexity pattern. If n is finite then = 00 otherwise Ni are all finite. The multiscale 
function of order m associated to A and N is the function defined by 



(9) 


where k satisfies Ni + - ■ •+iVfc_i < m and Ni + - ■ •+Afc > m and NJfi = m—{Ni + - ■ ■+Nk-i). 

We denote by Fa,n the limit of this function when m goes to infinity. Such a function is 
called a multiscale function. 

The three sequences A, N and T encodes the scale structure of a given multiscale function. 

As an example, we can define the Multiscale Okamoto’s functions as follows : 

Definition 8 (Multiscale Okamoto’s functions). — Let a = (ai,... , 0 ^) be a sequence 
(potentially infinite) ai g]0, 1[ and N = {Ni,... ,Nn) such that Ni € N*, i = 1,... ,n. If n is 
finite then Nn = 00 otherwise Ni are all finite. The multiscale Okamoto’s function of order 
m is the function defined by 



( 10 ) 


where k satisfies Ni + - ■ •+Afc_i < m and Ni + - ■ ■+Nk > m and NJfi = m — (Ai + • • •+Afc_i). 

We denote by Oa,N Ihe limit of this function when m goes to infinity. Such a function is 
called a multiscale Okamoto’s functions. 
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The three sequences a, N and T encodes the scale structure of a given multiscale Okamoto 
function. 

We denote by Okamoto the set of all multiscale Okamoto’s functions. In what follows, we 
abbreviate a multiscale Okamoto’s function by MSO function. 


2.6. Structure of MSO functions : A characterization using symbol seqnences. 

— In order to characterize multiscale Okamoto’s function with symbol sequences, we use 
some well known construction on space of symbol sequences as exposed in |321 §.2.2.a p. 96]. 

Definition 9 (Space of symbol sequences). — Let A he a subset of a given metric space. 
We denote by T,a the space of symbol sequences on A defined by an infinite Cartesian product 
of copies of A, i. e. 

OO 

^a = IIa. ( 11 ) 

i=l 

An element of T,a is called an infinite sequence or word defined on A and is denotes by 

s € TiA, s = (si,..., Sn, ■ • •) where Si € A, Vi € N. (12) 

Let a, N be two admissible sequences. We denote by Sa,N the infinite sequence in S[o,i] 
defined by 

■Sa.N = (ai . . . 0102 . . . 02 . . .) (13) 

where oi is repeated Ni times, S 2 is repeated N 2 times and so on. Using this construction, 
we have the following Lemma : 

Lemma 10. — The set Okamoto is in one to one correspondence with S[o,i]. 

2.6.1. Topology on O\isim.oto. — We introduce a topology on the set of multiscale Okamoto’s 
functions using the previous characterization. 

A topology on can be defined as follows. 

Definition 11 (Metric). — Let A he a metric space and d its metric. A metric on T,a is 
given for all s, s' € Tia by 

’ ^ ^2M + d(s„s')- 




We refer to [321 p.98]. 

The previous metric induces a topology on the set of multiscale Okamoto’s functions via 
the one-to-one correspondence with S[q 1 ] x Spj* as follows : 

Definition 12. — A metric on Okamoto is given for all /a,N o-n-d /a',N' € Okamoto by 

I^Okamoto(/a,N)/a',N') — ('Sa,N) ®a',N')’ (^h) 
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2.6.2. Symbolic dynamics and the three basic classes. — The study of multiscale Okamoto’s 

functions is related to the behavior of the curve under change of scales. Precisely, the data 
of the two sequences a and N encodes the way the curve behaves over the set of scales 
T = {TjjigN*- The change from scale T* to Tj+i is then associated on the sequences side to 
the classical shift map a : —>■ defined for all s G by |32l §.2.2.b p.lOO] : 

a{s)i = cTj+i, i € N*. (16) 

The shift map is the simplest example of a chaotic map (see [321 §.2.1.e p.93]). The main 
property of these maps is summarized in the following Theorem (see [321 Prop.2.2.11 p.l05]); 

Theorem 13. — The shift map possesses the following properties : 

1. a countable infinity of periodic orbits, 

2. an uncountable infinity of non periodic orbits, 

3. a dense orbit. 

This Theorem can be used to classify Multiscale Okamoto functions depending on their 
Scale structure. Precisely, we introduce the following classes of functions : 

- Self similar MSO functions denoted OkamotOseZ/ corresponding to periodic sequences. 

- Random MSO functions denoted by OkamotOj.a„c; corresponding to non periodic se¬ 
quences. 

- The Chaotic MSO function corresponding to the dense orbit. 

It must be noted that the scale structure of a given MSO function can be very complicated 
independently of its shape structure which is more related to its regularity. It means that the 
previous encoding does not capture the geometrical complexity of a MSO function. This is 
precisely the starting point of the scale dynamics studied in the following Sections. 

2.6.3. Complexity of Multiscale Okamoto’s functions. — There exists several notions of com¬ 
plexity related to symbol sequences in combinatoric. Using these notions, it is possible to select 
from a given symbol sequence a natural candidate with low complexity and representing the 
associated multiscale Okamoto’s function in one of the previous family. We refer to the report 
of J-P. Allouche in [2] for more details and precise definition of these notions. 

2.7. Possible generalizations ?— The previous definitions and constructions can of course 
be generalized in many directions. For example, the actual definition of multiscale functions 
use the same elementary action in order to produce the function between two successive 
scales. In order to define more complicated objects, one can allow a given set of elementary 
actions to go from one scale to the next one, producing mixed multiscale structures on a 
given function. The formalization of this kind of objects is of course a little complicated but 
does not bring new fundamental ideas despite its interest. As a consequence, we prefer to 
discuss the construction of our formalism in the limited (but already substantial) setting of 
multiscale functions. 


3. The modeling problem 

In the following, a multiscale Okamoto’s function is the mathematical analogue of a real 
physical process. It means that we assume that the exact behavior of the system is given by 
such a function. The modeling problem is then formulated as follows : 
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Modeling problem; Having a finite set of data concerning a given multiscale Okamoto’s 
function, can we recover the nature of the underlying funetion ? 

As we have already discussed in the first Section, the nature of an object is for example ; 
random or not ? or depending on scales ? etc. 

As we will see in the following, the answer is no. Only an infinite amount of data can give 
a complete characterization of the nature of an object. This result has strong consequences, 
as it means that no models can decide the exact nature of a physical phenomenon. In order 
to do so, we need some extra conditions which are coming from physics and more precisely 
from the underlying framework of a given theory. As an example in general relativity, the 
models are constructed assuming that no randomness occurs and that all the objects can be 
described in the differential framework. In other words, an identification of a given “exact” 
model is related to a given “philosophy of nature”. 

The previous modeling problem is not sufficiently precise in order to formulate a result. 
Using these results, we can precise the modeling problem for multiscale Okamoto’s functions. 

Modeling problem 2. Let m G N* 6e given and /a,N G Okamoto. Can we decide 
the “nature ” of /a,N knowing the associated multiseale Okamoto’s funetion of order m for 
arbitrary m ? 

3.1. Toward scale dynamics. — The previous discussion proves that the complexity of 
a given multiscale function is not related to the nature of the limit function. The basic 
information which is missing is the behavior of the seale derivatives (left and right) during 
the scale process, i.e. one needs to introduce a seale dynamies. This is precisely what is 
developed in the next Section. 


4. Scale calculus 

In this Section, we define the notion of seale derivative which will be used in the next 
Section in order to define what we call a scale dynamies. Our definition is based on some 
classical tools of discrete ealeulus (finite differences). We use the notations which are usual in 
the theory of the time-seale ealeulus indicating by this way that most of the notion defined in 
this section can be generalized over multiscale functions with a general time-scale sequence. 

4.1. Reminder about time-scale calculus. — We consider T a discrete and finite time- 
scale with a = min(T), b = max(T) and card(T) > 3. 

Definition I 4 . — The backward and forward jump operators p,a :T —)■ T are respectively 
defined for a// Vt G T by : 

p{t) = sup{s G T, s < t} and a{t) = inf{s G T, s > t}, 
where we put sup 0 = a and inf ^ = b. 

Definition 15. — The forward graininess (resp. backward graininess) funetion p : T —>■ 
(resp. V : T —)■ M+ ) is defined by p{t) = a{t) — t (resp. v{t) = t — p{t)) for any t gT. 
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We set = T\]/3(5),6], = T\[a, (T(a)[ and n T^. Let us recall the usual 

definitions of A- and V-differentiability. 


Definition 16. — A function 
differentiable at t G T'^ (resp. 


in 


u : T —M"' , where n € N*, is said to be A- 
S/-differentiable at t G T^) if the following limit exists 


lim 

s—yt 
S^(7{t) 


u{a(t)) — u{s) 
a{t) — s 


resp. 


lim 

S — 

s^p{t) 


u{s) - u{p{t)) 


s - p{t) 


In such a case, this limit is denoted by Au{t) (resp. Vu{t)). 


(17) 


Let us denote by f At the Cauchy A-integral defined in [Sj p.26] with the following result. 


Theorem 17 ( [Sj Theorem 1.74 p.27]). — For every u G there exist a unique A- 

antiderivative of u in sense of AU = u onT^ vanishing at t = a. In this case the A-integral 
is defined by 

U{t) = f u{t)At 
J a 

for every t gT. 


4.2. The scale derivative. — The aim of this Section is to define the natural object 
encoding the behavior of the A or V derivatives over the sequence of time-scales associated 
to a given scale function. 

Definition 18 (A Scale derivative). — Let Y be a given scale function over the scale se¬ 
quence T. The A scale derivative o/F is the scale function denoted by A(F) and is defined 
by : 

- = Tf. 

- [A(F)]^ = A(F,)GC(Tf,]R) . 

As an example, the scale derivative of the scale Okamoto function of order a = | is 
given by For multiscale functions, the scale derivative is a priori not a multiscale func¬ 
tion but only a scale function. This remarks put in evidence the disconnection between the 
complexity in scale of a given function and its complexity from the scale dynamical view point. 


A notion of V scale derivative is dehned in the same way. We denote by CaCIT) (resp. 
CvCT)) the set of functions which are scale A differentiable (resp. scale V differentiable) 
over T. 


4.3. The scale antiderivative. — Although we will not use scale antiderivative in a hrst 
approach to scale dynamics, we will need for some applications to integrate scale functions. 
As a consequence, we provide the corresponding notion. 


Definition 19 (A Scale antiderivative). — Let F be a given scale function over the scale 
sequence T. The A scale antiderivative of F is the scale function denoted by J F A and is 
defined over the scale sequence T by 




F,A€C(T„M). 


(18) 
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Figure 6. A Scale derivative 


As an example, the scale antiderivative of the scale Okamoto function of order a 
given by 




Figure 7. A Scale antiderivative 


4.4. Scale equations. — The scale derivative and antiderivative can be used to define 
what can be called scale equations which are more or less time-scale equations on each scales. 
These equations will have an important role in the last Section concerning applications of 
scale calculus and scale dynamics to classical mechanics. 

We define scale equations as follows : 
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Definition 20 (Scale equations). — A scale equation is a formal relation written on scale 
functions F on a given scale sequence T of the form 

RTfiFi,A, J ■A) = 0, (19) 

where F = {F)} and T = {Tj}. 

As an example, we can consider the linear scale equation given by 

A{Fi) - /i(Ti)F, = 0. (20) 

Some applications will lead to special scale equations which we call scale invariant and are 
defined as follows : 

Definition 21 (Scale invariance). — A scale equation is scale invariant if the scale equa¬ 
tion 

Rj^(^Fi,A, j -A^ =0, (21) 

where F = {F,} and T = {Tj} satisfies 

Rt, = R, ( 22 ) 

where R is a fixed relation. 

The previous linear scale equation (l20]l is not scale invariant. Indeed, the scale operator is 
given by 

Ft, (f, a, j = A(F) - (23) 

which is explicitly scale dependent through the graininess constant /i(Ti) which depends on 
each scale Tj. 


As an example of a scale invariant equation, we introduce the scale Euler-Lagrange 
equation which will be studied in the last Section. 


A Scale Euler-Lagrange equation is dehned for all F by 

/or \ or 

V ^(r..AF.) =-(F..AF.), 


(24) 


In this case the relation is given by the Euler operator using scale calculus and given by 

dL dL 

^Euler 

over scale functions whose form is independent of scales. 


(25) 


Of course, more complex notion of invariance can be defined. Eor example, in many physical 
problems we can expect only a partial scale invariance, meaning that the equation keep the 
same form only on a given range of scales. This problem will be discussed in details in Section 

E 


MULTISCALE FUNCTIONS, SCALE DYNAMICS AND APPLICATIONS 


17 


4.5. Scale equations and scale regimes. — The study of scale equations depends mainly 
as in the classical theory on the underlying functional space on which such an equation is 
studied. In the following, in particular in Section [71 we will deal with the dynamical behavior 
of scale equations over multiscale functions. In that case, the notion of scale regime will have 
a very important role. 


In order to identify scale regime for a given scale function, one need to select a given 
class of scale comparison functions. Examples are given by the classical Hardy scale., the 
logarithmic scale and the power law scale. We refer to [30i §.2.5 p.23] for more details. We 
only develop the power-law case in the following. 


The power-law comparison scale is defined by the family (see [301 p.24]) : 

V = {f^{t) = e, a>0}. (26) 

The exponent is obtained by looking for the quantity 

ln(/r) • ^ ^ 

In order to precise the scale regime, we introduce the definition of scale range: 


Definition 22 (Scale range). — Let T be a given scale sequence. A scale range between 
two time-scale T^o and T^i in T, with mo and mi G N such that mo < mi, is denoted by 
[TmojTTmi] and is defined by 

m\ 

[T™o,T™i]= U 

m>mo 


Let T be a given scale sequence and consider a scale function X over T. We then are lead 
to the following dehnition of a scale regime : 


Definition 23 (Pointwise scale regime). — Lett G T^o for mo > 0. The pointwise scale 
regime of X. in t, denoted by a(X,t), is the quantity defined by 


[a(X, t)]j^ 


ln(/r^|A(X)^(t)|) 

ln(/rm) 


(29) 


for all m > mo. 


Definition 24 (Local scale regime). — The local scale regime ofX, denoted by Q!(X), is 
defined by 

[a(X)]T„, = sup [a(X,t)]Tp , (30) 




for all m > mo. 


In applications, the important information is related to the evolution of the global scale 
regime when m goes to inhnity. Precisely, we have : 

Definition 25 (Global scale regime). — Let mo, mi G N such that mo < mi. The global 
scale regime ofX over [Tmo^Tmi] is defined by 

sup [a(X)]Tr^. 


( 31 ) 
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Remark 1. — These definitions are of course reminiscent of the classical pointwise and local 
Holder exponent for continuous functions as defined for example in |29] . Nevertheless, it must 
he noted that we do not need to assume that there exists a continuous function associated to 
the given scale function in order to defined such exponents. They are constructed directly on 
the family of discrete data which is, from our point of view, the only information that one can 
obtain in experimental settings. 

Example : Consider the multiscale Okamoto’s function defines by a = {2/9, 2/3, 5/6} and 
N = {4,3,oo}. Its pointwise scale regime over the scale range [Ti,Tio] is illustrated in the 
Figure [51 It is computed for the point to and to,i in Ti as by construction, the pointwise scale 
regime of to ,2 coincide with the one of to- 



(a) /imAOa,N with respect to in logarith¬ 
mic axis. 


aTS;i(Oa,N,to) 


Ti T2 T3 T4 Ts Te T7 Tg Tg Tjo 

(b) Pointwise scale regime 


Figure 8. MSO function with a = {2/9, 2/3, 5/6} and N = {4, 3, 00 } 


On this simple example, we can see the three pointwise scale power law regime over the 
scale range [Ti,Tio]. We also display the power law function with ai = 02 = 

and 03 = ^00^ in order to compare to the slope of the pointwise scale regime. In that case, 
the pointwise scale regime is the same for the two points. In consequence, the global scale 
regime is a power law regime with a = a^. 


In what follows, for a general scale regime, we denote by (X) the global scale regime 
of X (or simply scale regime) over the scale range [TmojTmi]- 


5. Scale dynamics 

A scale function being given we are interested in the following problems : 

- Can we determine if the scale function is a multiscale function ? In this case, try to 
determine the elementary scale actions. 
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- Assuming that we know a given scale function up to scale m G N*, can we precise the 
“best” in some sense continuous limit model ? 

To do that, we see that it is necessary to understand the behavior of the scale derivative 
under change of scale, i.e. scale dynamics. 

We begin first by studying multiscale Okamoto’s function and we prove that the scale 
pattern complexity can be recovered by looking at the dynamics of the scale derivative. 
Then, we provide general transformation formula for the scale derivative between two given 
time-scale. 

5.1. Scale dynamical analysis of Multiscale Okamoto’s functions. — In this section, 
we study the behavior of the A-derivative of multiscale Okamoto’s functions under change of 
scale. 

5.1.1. Scale dynamics of Okamoto’s functions. — We begin with single scale Okamoto’s 

function. Considering Tq = and a €]0,1[. By definition of the scale Okamoto’s 

action Oa (see Definition [5|) , we obtain the time-scale Ti = To,Oa = {^Oj ^ 0 , 1 ; ^ 0 , 2 , ii} where 

to,i = to + for f = 1, 2 and = ti — to. Considering the discrete function Oa,To defined 

O 

on To and the discrete function Oa,Ti defined on Ti with the points obtained by the scale 
Okamoto’s action. 


We introduce the notion of reference scale regime which we choose to be the linear scale 
regime, meaning that the dependence of the correction term with respect to Hi is linear. This 
assumption leads to the notion of discrete reference function on Ti associated to a given 
discrete function on Tq, denoted by [OoIto- 

Remark 2. — The choice of the linear scale regime is of course arbitrary from the math¬ 
ematical point of view. However, from the physical side, it corresponds to the fact that we 
always compare a possible smooth one which is the standard setting in Physics. Any deviation 
to this framework is interpreted as the need for a new mathematical setting (non-differentiable, 
stochastic, etc). 


iTi 


For Okamoto’s function, the discrete reference function on Ti is defined by [OaJ-fo lTo = 
Oa,To and 

Oa,To(^l) — Oa,To(^o) 


[^a]To(fo.l) = Oa,To{lo) + Ml 


Mo 


[^a]To(fo.2) — Oa,To(^o) + 2M1 
By definition, we have 


Oa,To(^l) — Oa,To(^o) 


Mo 


AOa,To(to) = AOa,Ti(to) = —AOa,To(io) 

tl — to p. 


then, we obtain : 


Proposition 26. — The A-derivative of Okamoto’s function satisfies 

AOa,Ti(io) = A[Oa]^J(to) + C(to) 
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where 

C(io) = ^OaJoito), 

In order to generalize our previous results, we introduce the following definitions : 

Definition 27. — Consider a scale function F on T. We denote by F* the reference scale 
function associated to F defined at scale by ^ for all m> 1. 

Definition 28. — Let A be a given scale operator on scale functions. We call scale effect 
induced on A the difference between the action of A on a given scale function and its 
associated reference scale function. 

Remark 3. — The quantity C{to) in Proposition [2^ corresponds to the scale effect on A 
between Tq and Ti. 

Consider now the scale sequence T associated with the scale or multiscale Okamoto’s func¬ 
tions. We have ; 

Proposition 29. — The scale effect induced on the scale /S.- derivative of Okamoto’s function 
is given by 

AOa = AOl + C (Oa) 

where 

C (O.) = (^ -1) Ao;, 

Analogous formulas holds for the scale V-derivative. 

From this proposition, we recover the features of Okamoto’s function given in the Theorem 
[1] assuming p = 1. 

5.1.2. A scale dynamical approach to the Okamoto-Kobayashi’s theorem. — The analytic 
properties of the Okamoto’s function are related to the asymptotic behavior of the A and V 
derivatives and the correction term with respect to scale. Indeed, the existence of a derivative 
for a given point t G [0,1] can be check as follows : assume that t € for a given m G N. 
The limit function admit a derivative if and only if we have 

- lim [A(Oa)]m,(t) and lim [V(Oa)]m(i) exist. 

m—)-oo m—)-oo 

- lim [C{Oa)]m = 0 . 

m—)-oo 

The hrst condition implies that the left and right derivatives exist at point t for the limit 
function and the second condition implies the equality of the left and right derivative which 
implies derivability at point t. 

An easy case is obtained when the previous conditions are satisfied from a given scale. In 
particular, we have : 

Proposition 30. — The correction term C (Oa) vanishes if and only i/a = 1/3 or a = 1/2. 

Proof. — The first point follows directly from equation ()29p in Proposition 1291 

The second case follows from the following observation : in order for a given scale to 
have a point such that AF = 0, one need to construct in the iterative procedure for a given 
elementary time scale as given in Definition [5] a conhguration such that Oa (io,i) = Oa(to, 2 )- 
This is possible if and only if a = 1/2. □ 
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The situation between the case a = 1/3 and o = 1/2 is nevertheless very different. 

- For 0 = 1/3 the correction term vanish in all point of the scale sequence. This result 
induces the third point of the Okamoto-Kobayashi’s theorem. 

- For a = 1/2, the correction term vanish only on the following left reduced scale sequence 

obtained from the Okamoto’s scale sequence taking at each step from the action on an 
elementary time-scale Teiem = the point to,i from Teiem,Oa- We denote by 

this scale sequence. 



Figure 9. The points in red correspond to the image of the points by Okamoto’s 
function 

We have for all t G there exists m{t) G N such that t G and [A{Oa)]mit) = 0 for 
all m > m{t) . 

This result extends easily to the V case changing the left reduced scale sequence for the 
right scale sequence obtained in the same way taking to ,2 instead of to,i in the iterative 
construction. 

It must be noted that we have no information on the point of the sequence belonging to 
the Cantor’s set T \ U and defined by 

oo 3P-1-1 

Cantor = [^ 

p=l k=0 

In order to conclude about the derivability of the limit function, one needs to check that 
the A and V derivatives possess a limit when m goes to infinity. We will not reproduce that 
computations made for example by H. Okamoto or K. Kobayashi in their papers as this is 
not the purpose of our work, but one can of course interpret their results in our setting. 

5.1.5. Finite scale sequence and asymptotic model identification. — Let us begin with a 
definition : 


3A: + 1 
3P 


U 


3k+ 2 
3P ’ 
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Figure 10. Correction term for Okamoto’s function a = 1/2 


Definition 31. — LetY he a scale sequence and F a scale function. We call F^o an asymp¬ 
totic model for F. 

For example, Oa is an asymptotic model for Okamoto. 

Identification problem : Assume that we have access to F up to the time-scale T^- 
Can we characterize the asymptotic model with these data ? 

In the Okamoto case, one needs only to identify the parameter a. This parameter is 
completely fixed using a finite set of data. As a consequence, the asymptotic model is 
determined knowing only a partial set of observations. Of course, this result is due to the 
fact that we presuppose that the scale regime will not change, i.e. that o is independent of 
scale. This assumption can be though as an a priori point of view on the nature of the real 
phenomenon. 

This result will drastically change in the multiscale case. 

5 . 1 . 4 . Scale dynamics of multiscale Okamoto’s functions. — The identification problem can 
be more complicated than the one which is exhibited by the Okamoto’s function. In the 
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multiscale version, the problem is worse. 

We can perform the same kind of computations for a multiscale Okamoto’s function Oa,N 
with a = (ai,..., an) and N = {Ni, ..., Nn) as in Definition [8l Indeed, following the same 
idea as previous with the scale Okamoto’s functions, we obtain the following proposition : 

Proposition 32. — The scale effect induced on the scale A-derivative o/Oa.N is given by 

AOa,N = AOa,N + C (Oa,N) 

where 

C(Oa,N)= AO*n 

and where Sa,N is the sequence defined in Equation Ii31 An analogous formula holds for the 
scale V-derivative. 

Here again, the asymptotic properties of the model are fixed by the set a and N. However, 
the identification problem can not be solved. 

Indeed, assuming that we know the function F up to scale m, we do not know, if the 
sequence Sa,N is hnite or infinite. Even in the finite case, we do not know the length of the 
sequence. Moreover, even if we know this length, let say k, and we identify already k — 1 
terms in Sa,N up to scale m, we do not know when the next change will produce. 

This means that, without an assumption on the asymptotic model, this is impossible to 
conclude from a finite scale sequence observation. This assumption is not induced by the set of 
accessible data and can not be checked. It comes necessarily outside of the given framework. 

Remark 4- — In more physical terms, the choice for a particular model can then be only 
justified using a certain philosophical point of view on the nature of a given phenomenon which 
can not by definition and construction be proved by any experimental devicel^. 

One must noted that the previous remark is in fact valid for all kind of model which can 
be constructed from a given set of experimental data as, in practice, this is impossible to have 
access to the full scale seauencew'^\ A limit model can not and will never represent the reality 
of a phenomenon as we have no possibilities to select between different admissible models up 
to a given scale. 

5.2. General formulas for the scale effect on A and V. — 

5.2.1. Passage from a scale Tq to a scale Ti. — In this section we provide a general 
transformation formula for the scale derivatives between two given time-scale. 

We want to define relations between two discrete and finite time-scales Tq and Ti in order 
to study the effect of changing scale on the scale A-derivative and the scale V-derivative for 
a given scale function. We are interested in the case where Tq C Ti with Tij-j-o = Tq and 
we consider the case where Tq and Ti are uniform which means their graininess function 
are constant. Moreover, without loss of the general idea, we consider the case where the 
graininess function pQ of the time-scale Tq and the graininess function pi of the time-scale Ti 


1. This is the case for example of String theory which assumes that physical properties are coming from 
the description of a particular geometric object generalizing Einstein’s point of view on space-time. 

2. In Physics, we have for example the well known limitation due to Heisenberg uncertainty principle. 
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are such that fiQ = 2/ii. With such considerations, we put Tq = {toj^i} and Ti = 
We consider the discrete reference function defined by 


FS(t*) = FTo(to) + 


To 


-^To(^l) - -^To(io) 


^^0 


Ail- 


If we compare (t*) with F^^ (A*), we can write AFTri(Ao) as 


AFtiC^o) = 


Ft, (a*) - Fj (A*) + F^^i (A*) - Ft, (Aq) 


To 


Ail 


We obtain, using the definition of F.^^^ (A*) 


FT,(A*)-F;i(i*) 


AFt,(Ao) = AF;^i(Ao) + ^ 


Ail 


(32) 


(33) 


(34) 



Figure 11. Comparison between the expected point and the new point 


The explicit dependence of the scale behavior is shown with the last term on the right. In 
order to express this term as a dynamical quantity, that is to say in mean of A and V, we 
note that by definition of the A-derivative, V-derivative and F.^^^(A*), we have 


Ft,(A*)-F5(A*) (V-A)^ 

(35) 

Ail 2 

Finally, we obtain 


AFt,(Ao) = AFrj,^ (Ao) + 2 -^Ti o Ti(io)- 

(36) 

Doing the same with the V-derivative, we obtain 


VFt,(Ai) = VFJ((Ai) + (^^Ft, opi(Ai). 

(37) 
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As we see, the deviation from the reference function at scale Ti is governed by the difference 
between the A and V derivative of the -Fti • We give a general version of these computations 
in the next Section. 

5.2.2. Scale dynamics. — Consider now the scale sequence T where each time-scale Tj G T, 
i > 0 , has its graininess function fii such that /r* = 2 /Xj+i where is the graininess function 
of the time-scale Tj+i. 


As we have seen in our previous computations, an important role is played by the operator 
V - A 

—-—. In particular, the correction term can be expressed as an action of this operator. In 
order to simplify the formula in the general case, we introduce the following notation : 


Definition 33 (Correction term). 

V - A 

tion terms defined by Cx (F) = —-— 


We denote by and Cy, the left and right correc- 
V-A, 


x(F)] andCy{Y) = 


[xi(F)] with 


and 


[x(F)1t„ = 




Fop 

over 

1 

n 

(38) 

F 

over 

Tm \ Tm-1, 

F o (T 

over 

Tm -1 c 

(39) 

F 

over 

Tm \ Tm-1, 


For convenience, we introduce the following scale sign function : 
Definition 34- — We denote by e the scale sign function defined by 


— 

'■ -I h 771, 


-bl over Tm-i C T^, 
-1 over Tm\'irm-i. 


(40) 


Considering the scale reference function F*, we have the following proposition : 

The scale effect induced on 


Proposition 35 (Scale effect on the scale derivatives). 

the scale A and V derivatives of F is given by 

AF = AF*+£-Cx(F), VF = AF*+£-Cx(F). 


(41) 


In many applications, we will need to have a convenient formula for the scale effect induced 
on the A and V derivative acting on a scale functional, i.e. a functional defined on a scale 
function. This is provided by the following Lemma : 

Lemma 36 (Scale effect and chain rule). — Let / :MxM— 7 >M 6 ea sufficiently smooth 
real valued function. We have for all scale function X on the scale sequence T, the following 
formula 


j-i 


A/(T,X) = A/(T,X*) +e • ^ ^[Cx(X)p|^(x(T), x(X*)). 

i>i 


(42) 


and 


V/(T,X) = V/(T,X*) + e • ^ Li)L^[C,(X)p|^(x(T), x(X’^)). 


i>i 


dx^ 


(43) 
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Proof. — We detail the proof only for A as the computations are equivalent in the V case. 
We fix a given scale and we need to distinguish between points in T^-i C and those 
in Tm \ Tm-l- 


For t G Tm\Trm-i, we have X.oa{t) = 'K*oa{t) by construction, so that f(Toa(t),'Koa{t)) = 
/(T o o a{t)) and X(t) = X* o a{t) — ^AX(t). As a consequence, we obtain 

^iAf{T, X)(t) =/(T o a{t),X* o a{t)) - /(T(t), X(t)), 

=/(T o a{t), X* o a{t)) - /(T(t), XV(t) - ^AX(t)), 

=/(T o a{t), X* o a(t)) - /(T(t), XV(t) - ^AX’^(t) - (X)(t)), 

=/(T o a{t), X* o ait)) - /(T(t), X*(t) - (X)(t)). 


For t G Tm \ Tm-ii we have £(t) = — 1 so that 

mA/(T, X)(t) = /(T o C7(t), X*(t)) - /(T(t), X*(t) + (X)(t)). (44) 

As / is sufficiently smooth, we can make a Taylor expansion with respect to (T(t), X*(t)). 
We obtain 


mA/(T, X)(t) =/(T o a(t), X* o a(t)) -/(T(t), A*(t)) - ^ (X)(t)F|^(T(t), X*(t)). 

i>i ^ 

(45) 

We then obtain for t G \ Tm-i 




A/(T,X)(f) = A/(T,X*)(f) +eit) • ^ ^[C,(X)(t)]^|^(x(T)(f), x(X*)(t)). (46) 

i>i J- 

For t G Tm-i C Tm, we have X(t) = X*(t) and X o cj(f) = X*(f) + fiA'Kit). As a 
consequence, we obtain 


/iA/(T, X)(t) =/(T o a(t), X o c7(t)) - /(T(t), X*(t)), 

=/(T o a(t), X*(t) + /iAX(t)) - /(T(t), X*(t)), 

=/(T o ait), X*it) + /iAX*(t) + ^£(t)C>, (t)) - /(T(t), X*(t)), 
=/(T o ait),X* o a(t) + fieit)C^ it)) - /(T(t), X*(t)). 


For t G Tm-i C Tm, we have e(t) = 1 so that 

/iA/(T, X)(t) = /(T o ait), X* o a(t) + /iC,(t)) - /(T(t), X*(t)). (47) 

As / is sufficiently smooth, we can make a Taylor expansion with respect to (Tocj(t), X*ocj(t)). 
We obtain 


MA/(T,X)(t) = /(Toa(t),X*oa(t))-/(T(t),X*(t))+^^[C>,(X)(t)P^(Toa(t),X*oa(t)). 

i>i * 

(48) 


We then obtain for t G T^-i C 




:9V, 


A/(T,X)(t) = A/(T,X*)(t) +£(t) . ^^[CVX)(t)]^'^(x(T)(f), x(X*)(f)). (49) 

i>i 


□ 



MULTISCALE FUNCTIONS, SCALE DYNAMICS AND APPLICATIONS 


27 


6. Asymptotic differential operator 

The previous scale analysis can be used to define a new class of differential operator on 
the asymptotic class of functional set covered by limit of scale functions. These operators 
generalize the Box derivative introduced in m et the Ito calculus for stochastic processes 
(see [27j ) . The basic idea is to used the natural decomposition of a given scale function having 
a fixed scale regime in its “regular” part and its deviation part which is responsible for the 
scale regime. 

6.1. Extension and decomposition of scale functions. — Let mo < mi and [TmoiT^i] 
be a scale range over which a scale function X possesses a scale regime of order 0 < a < 1. 
Then, we have : 

- We call extension of X and we denote by ^ext a scale function such that X | ^ 

and X-ext possesses a scale regime of order a over [T^i+i, Too]- Then Xext,oo belongs to 

- Moreover, we have of order for all m > mo and j > 1. 

- By construction, we have [X*]ea;t which possesses a linear scale regime and is associated 
to a piecewise function denoted by X^. 

As a consequence, we can write X^xt = ^txt^^ext where De^t = Xext—X*^^. The asymptotic 
limit of X*^^ denoted X^o is decomposed as 

Xoo = X^ + Dao, (50) 

where X^ is a piecewise differentiable function and the deviation part from this differentiable 
behavior is given by D^o which is a particular way of decomposing X^o into a regular and 
non regular part 

Remark 5. — In [10] this operation is made directly on a very special functional space for 
which the non regular part is fixed in a given class. This allows the authors in this case to 
defined a projection which gives well defined regular and non regular part. 

6.2. Asymptotic differential operator. — The A and V derivatives are asymptotically 
equivalent to the classical derivative as long as the scale regime is linear. However, they do 
not possess asymptotic limit on scale function possessing a scale regime of order 0 < a < 1. 
Using the previous extension/decomposition of a scale function over a given scale regime, we 
give a meaning to a differential operator acting on Aqo by taking advantage from the fact 
that a part of the decomposition admits a classical derivative. 

Definition 37. — Let X be a given scale function and RTrnQjm^ ® scale regime with mo < 
mi. We denote by X^xt, X*^^ and D^xt the extension and decomposition of X associated to 
RTrnoJmi X^o, Xf^, Doo thc osymptoUc limit and decomposition ofX. We then define 
the operator Aqo acting on Xoo as 

Aoo[Aoo] = (51) 

As a consequence, the operator Aqo extracts the derivative of the regular part of Xoo. 
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Another way to formalize this operation is to introduce the operator reg on scale function 
such that reg(Xea;t) = In this case, we have 

Aoo[Xoo] = Aoo[reg(Xe,.t)] = (52) 

Remark 6. — This definition must be compared with the definition of the Box derivative in 
[IQ] and also the definition of the Nelson’s forward and backward derivatives over stochastic 
processes in |24j . 

For applications, we need as usual the behavior of this operator over composition of func¬ 
tions. In order to have explicit formula we introduce a special functional space over scale 
functions : 


Definition 38. — Let X be a given scale function and a power law scale regime 

of order 0 < a < 1 with mo < mi. We denote by X^xt, ^ext the extension and 

decomposition of X associated to We denote by C'a,a+ (resp. the set of 

scale functions such that for j a = E{l/a) there exists A+ (resp. A+j such that 

lim /i^^-^[C>,(T,Xe.t)p = A^ (resp. lim X^xtW^ = xt) . (53) 

We have the following result : 


Proposition 39. — Let X be a given scale function and R^^ ® power law scale regime 

of order 0 < a < 1 with mo < mi. We denote by X^xt, D^xt the extension 

and decomposition ofX associated to and Xoo, Xf^, the asymptotic limit and 

decomposition of X. Assume that X € C'a,a+- We then define the operator A^o acting on 
f{t,Xoo) as 


r/+ A-^,“ BN f 

A^fifiX^) ■.= ^fit,x:^) + 


(54) 


Proof. — We have 


Af{T,X,xt) = A/(T,X:,,) +£ • ^ ^[C,(Xe.t)Pp^(x(T), x(X:,,)). (55) 






i>i 


dx^ 


j-i d^f 

As the quantity [Cxi(Xea;t)]-^ j(x(T), ^(Xg^,^) admits a linear scale regime for j > ja 
by assumption, we have that 


reg 




^[C,(X„,)pA((x,(T), N(x:,,)) 

i>i 


ji ' ' ' dx^ 




^;^[c,(x„,)F^(>.{T). >.(x:,,)). 

j>jc 

(56) 

However, the asymptotic of these quantities are trivial for j > ja so that the limit reduces to 

1 ■ OE f 

-^Ai“^(t,A^). 


ja'. dxN 


( 57 ) 

□ 


In the same way we have 
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Proposition 40 . — Let be a given scale function and ® power law scale regime 

of order 0 < a < 1 with mo < mi. We denote by LLext, '^txt L)ext the extension 
and decomposition of X. associated to o,nd X^o, Xf^, D^o the asymptotic limit and 

decomposition of X. Assume that X G Cs/,\_- We then define the operator V^o acting on 
f{t,Xoo) as 

V^f{t,X^) := -^fit,X*^) - (58) 


6.3. Extension to complex scale functions and comparison with the Box deriva¬ 
tive. — Following the strategy exposed in [6] and further developed in m. we can extend 
the asymptotic scale operator to complex scale functions assuming linearity of the operator : 

Definition 4^ (Asymptotic Box derivative). — The asymptotic Box derivative, de¬ 
noted by Doo is the linear operator defined over complex valued scale functions by 

1^00 = — (Aoo + Vcxd) + (Aoo — Vcxd) , (59) 

where = —1 and g = {—1,1, —i, i}. 

In [10] , the definition is very similar and is in some sense equivalent even if not formulated 
in the same formalism. The main problem in HQl was to extract from the one parameter 
family of averaging of a given continuous functions, some information which can be encoded 
in a kind of derivative. The one parameter family associated to a given function can be 
clearly seen as a kind of scale function. The scale regime is then responsible for the divergence 
of the left and right derivative which is taking into account in m using some projection 
operator on the set of convergent function depending on the parameter. This is the role of 
the regularization operator in our definition. 


Nevertheless, we point out that the main point in our approach is that no asymptotic 
object need to exist ! Te asymptotic object is constructed using the extension procedure and 
does not necessarily corresponds to the real asymptotic (if any) of the scale function. 

Using Proposition (391 and liOl we have : 


Proposition 42 - — Let X be a given scale function and cl scale regime with mo < 

mi. We denote by X^xt, ^ext the extension and decomposition of X associated 

to ond Xoo, X^, Doo the asymptotic limit and decomposition of X. Assume that 

X € C'a,a+ FiCv,a_- We then define the operator Deo acting on f{t,X^) as 


□oo/(i,Aoo) 


□ 


f{t,Xo^) + 


Xg d^-f 
jal 


ifiX, 


* ] 
00 J 


where 


+ irj (^A^“ + Ai“) 


(60) 


( 61 ) 


One can compare this result with the corresponding result in m- 
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7. Asymptotic models of scale equations 

The previous formalism can be used to derived what is the adapted asymptotic model 
associated to a given scale differential equation. 

Definition 43. — We call scale differential equation any formal relation of the form 

F(T,X,A,V) = 0. (62) 

As an example, a particular class of second order scale differential equation is given by 

VoAX = /(T,X). (63) 

These models are usually constructed from a given macroscopic model corresponding 
for us to an initial scale Tq. This macroscopic model is in many cases associated to a classical 
differential or partial differential equation like for example the classical Newton’s equation in 
classical mechanics. This representation must be thought as a particular asymptotic model 
assuming a very special scale regime, namely the linear one. 

We give two examples which will be used in the Section concerning applications of this 
idea. 

7.1. Linear scale regime. — Take a scale sequence T and the set of scale functions having 
a linear scale regime. In that case, using the scale dynamics equation, we observe that 

lim Cx(X)= lim Ck(X) = 0, (64) 

m—>-oo m—>-OD 

When m goes to infinity, we have by construction that goes to a closed interval [a, b] and 
any discrete function belonging to the functional set (^(Tm,^) is converging to ^([a, 6],M). 
The operator A and V are converging to the classical right and left derivatives. 

For the previous scale equation, we must assume that X G Cl, ^(T), i.e. X is scale 
A differentiable over T and AX is scale V differentiable over T. If X satisfies moreover 
that limrn^oo A(X)'f^ and limm->-oo V(X)'f^ exist, then the condition of linear scale regime 
implies that X is converging to Xqo G (^^([a, 6], M). 

As a consequence, under the linear scale regime assumption, the asymptotic model asso¬ 
ciated to the previous second order scale differential equation is simply the classical second 
order differential equation 

^ = /(^>a;). (65) 

Of course, this assumption is known to be false. The scale regime, even in the classical setting, 
changes drastically when going to the microscopic scale. The question is then to understand 
what is the correct asymptotic model when the scale regime is changing. Our formalism 
provide such an answer, taking always the deviation with respect to the linear scale regime as 
reference. In other words, this means that we compare the behavior in different scale regime 
with respect to the macroscopic one took as a reference. 
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7.2. Fractional scale regime. — If we replace the linear regime by a fractional one of 
order 0 < a < 1 the situation changes drastically. 


If the initial scale equation is equivalent under the change of variable / to the equation 

A/(T,Xe.t) = G(T,X:,,). (66) 

Then its asymptotic continuous model is given by 

Aoo/(t,Xoo) = G'(t,X^), (67) 

which can be written as a classical (partial) differential equation 


1 ■ f 

^/(*. vjo) + = 0(<. 

We can resume the previous approach by the following diagram : 

A/(T,X) = A/(T,X*) + C7e(T,X) 



^ext) 


X*^) = Xi) + Coo,7^(^, x^) 


( 68 ) 


Figure 12. Diagram showing the approach of the extension and the passage to an 
asymptotic model 

The important point is that equation (I68p is a classical (partial) differential equation. As 
a consequence, this construction allows us to interpret perturbations of classical models as 
effects due to different scale regime over a single scale invariant equation. An example of this 
situation is given in the next Section with the Newton’s equation. 

7.3. Asymptotic versus scale equations. — The natural object coming from modeling 
in Physics, Biology or any modeling based on experimental data can be formalized in the 
framework of scale equations, dealing only with a finite number of quantities. However, most 
of the present models are formulated in the context of the differential calculus or different 
generalization like the stochastic calculus of Schwartz distributions. As we have seen, these 
continuous models corresponds in our language to special asymptotic model of a given scale 
system. The asymptotic character of these models is in general completely lost in the current 
formulation of Physics, etc, which can be called the fundamental continnous modeling 
assumption which is, by no way, explicit in the existing literature. 
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Fundamental continuous modeling assumption : Assume that a given system exhibit 
a given scale regime. A continuous model is obtained assuming that the observed seale regime 
is in faet valid up to infinity. 

This fundamental assumption is very useful when one is trying to interpret some very 
strange result coming from the analysis or numerical analysis of some models. This is in 
particular clear in the context of fluid mechanics where the governing equation are precisely 
written under an asymptotic assumption which relate the macroscopic and microscopic 
behavior. 

Assume, to see the difficulties induced by this assumption, that the scale regime change at 
very small scale. Then, doing numerical simulation of the corresponding continuous model, 
one must wait for the emergence of nonphysical or strange solutions. The same is true for 
the analysis of current continuous models, when one studies solutions exhibiting some scale 
structures. The Navier-Stokes equation or Euler equation are in this respect very illustrative. 
Simulations of the Navier-Stokes equation induce most of the time unwanted solution with 
no physical meaning. The problem is then to choose between the validity of the equation 
itself and the validity of the numerical method. Some paradoxical solutions are known for 
the Euler equation which have no physical meaning but can be explicitly constructed. Eor an 
overview of results and questions in fluid mechanics, we refer to [1611^. The construction 
makes use of a recursive and scale behavior. The question is then to choose between again 
the equation or the definition of generalized solutions that one must introduce to cover for 
example turbulent behavior. 

Again, all these problems are due to the fact that an asymptotic model has been used base 
on the fundamental continuous modeling assumption which is not explicit in the presentation 
of these models. A given continuous model is valid as long as the fundamental continuous 
modeling assumption can be in some extent justified. This is not the case for fluid mechanics 
for example, at least for turbulent regime where scale structures are assumed to appear at all 
scales. 


8. Applications to partial differential eqnations : the Diffnsion equation and the 

Schrddinger equation 

If different scale regime are present, we obtain with the previous formalism different asymp¬ 
totic models with different correction terms view as deviation from the linear scale. This is 
illustrated as follows: 

Such a result is particularly interesting when studying a scale invariant equation. In¬ 
deed, a single scale invariant equation will provide a possibly infinite number of asymp¬ 
totic continuous models corresponding to different scale regimes. The fact that different 
continuous models exist for a given physical phenomenon is then understood as the fact that 
the phenomenon exhibit different scale regimes during scaling. As a consequence, the multi¬ 
plicity of continuous models does not imply that no universal equation underlies the 
phenomenon in the scale framework. Each continuous equation is valid in its own domain 
of scales. 
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^(X>fR.Q — 


+ 

-^oo,7?.o 

^oofRx — 

Xt, 

+ 


^oo,Tlk — 

Xt„ 

+ 

DoofUk 


Figure 13. Illustrations of different continuous models with different scale regime. 

The part in blue corresponds to the classical or regular part. The part in red cor¬ 
responds to all the deviations to the regular part depending on the scale regime TZk 
chosen. 

8.1. Scale Newton’s equation and the diffusion equation. — We consider the classical 
equation obtained by Newton to describe the dynamical behavior of a particle of mass m under 
the action of a force deriving from a potential U. Precisely, we call Newton’s equation the 
following ordinary differential equation of order 2 

^ = U'{x). (69) 

The scale Newton’s equation is defined by 

VoAX = 17'(X). (70) 

The choice of this equation as a scale analogue of the classical Newton’s equation is sup¬ 
ported by the following results : For each scale, this equation corresponds to the variational 
embedding in the sense of [12j of the classical Newton’s equation which means that the 
solution of the scale Newton’s equation coincide with the extremal of the scale embedding of 
the classical Lagrangian structure under the time-scale calculus of variations. 


Remark 7. — The class of scale equations dSl corresponds to the choice of A as the deriva¬ 
tive of the position. As the scale Newton’s equation dZQ]) is defined in this class, we have also 
its V version given by 

AoVX = t/'(X). (71) 


Let us assume that there exists a function 'tp{t,x) such that 


and 


AX = -27 


dx 


(T,X), 


AX = -VX. 


(72) 


We then have the following theorem : 


(73) 
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Theorem 44- — Assume that X possess a scale regime of order \l‘l. Then the asymptotic 
continuous model of the scale Newton’s equation under the change of variables if is given by 


dip 



1 

ip 



Xl d^ip 

2 dx'^ 


+ i^Uip = 0 . 

27 


(74) 


Specializing to the subspace defined by 


7 = - 



(75) 


we deduce an interpretation of the diffusion equation as the dynamical behavior of the classical 
Newton equation on a particular scale regime : 


Theorem 45 (Diffusion equation versus Newton equation) 

Assume that X possess a scale regime of order 1/2. Then the asymptotic continuous model 
of the scale Newton’s equation under the change of variables ip is given by 


dip Xf d'^ip , 1 , 


(76) 


Proof — By Proposition [321 we have 




d-xf.it) d fdHip)\ 


+ 




\ y 


it,xut)) 


dt dx \ dx 
1.2^/ dln{iP) \ 

-5x2 V 


As AX = —VX, we obtain 


d-Xp^{t) 

dt 




(t) ^_dln{iP) 


(i)). 












MULTISCALE FUNCTIONS, SCALE DYNAMICS AND APPLICATIONS 


35 


Replacing in the equation, we deduce 




-R: 


dx 

9^ / 51n(V’)\ 
dx"^ \ y 




d ( ain(V’)\ , * , ,, _ d 




1 f dip'V 
y9xy 




A 

dx 




dt 


\ dx 




As a consequence, developing the derivatives of in “ip, we have 


-A 

dx 


1 dip 
Ip dt 


1 


dip 




1 d^p 


^p‘^ \dx 

A? 


1 


1 d'^ip 


1 


^p dx'^ ip‘^ \ dx 


dip 


Ip dt ^ \ ^ 2 j Tp'^ \dx j 2 Ip dx'^ 


dip 


A? 1 d'^ip 


{t,xp,{t)). 


{t,x:,it)). 


Using the scale Euler-Lagrange equation, we deduce that 


A 

dx 


1 dip 
Ip dt 


2 


7 + ~ ) :^ ( ~ ) “ 


dip 

dx 


A?_ 1 d^ip 1 

-— H- U 

2 Ip dx"^ 27 


(t,A^(t)) = 0. 


As a consequence, we obtain the following partial differential equation 

^2 \ 1 /aj.\ 2 ^2 o2. 


dip 

This concludes the proof. 


Ai\ 1 /aV’\ Xtd^iP 1 




+ 2^^'*’ ^ “■ 


(77) 

□ 


8.2. Box scale Newton’s equation and the Schrodinger equation. — In |10] . the 

Box derivative allows to recover the nonlinear and linear Schrodinger equation. This is also 
the case here using the definition of the Box derivative in the context of scale dynamics. As 
in the previous section for the diffusion equation, let us assume that 


AX = -AX. 

and there exists a function ip{t,x) such that 

AX = - 27 ^(T,X), 


(78) 

(79) 
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Using the definition of the Box derivative, it leads to 


at 


dln('ll;) 


it))- 


(80) 


Combining the A and V version of the scale Newton’s equation with the assumption given 
in the Equation (I78p . we obtain what we called the Box scale Newton’s equation as 


□ 

□t 



U'(X). 


(81) 


Theorem 46. — Assume that X possess a scale regime of order 1/2. Then the asymptotic 
continuous model of the Box scale Newton’s equation under the change of variables if is given 
by 


— 2^7 




1 f dipX 2 
ijj \dx ) 2 dx'^ j 


+ Up; = 0 . 


(82) 


Specializing to the subspace defined by 

7 = ^, A2_ = \l = h^ (83) 

and with 7 = — 1, we deduce an interpretation of the Schrodinger equation as the dynamical 
behavior of the classical Newton equation on a particular scale regime : 


Theorem 4 7 (Schrodinger equation versus Newton equation) 

Assume that X possess a scale regime of order 1/2. Then the asymptotic continuous model 
of the scale Newton’s equation under the change of variables is given by 


df) h^ 


Up). 


(84) 


Proof. — The proof follows from the same kind of computations as for the diffusion equation 
only replacing Aqo by the asymptotic Box derivative. □ 


9. Conclusion 

From the modeling problem, we develop in this paper a formalism called scale dynamics, 
which allows us to deal with multiple finite and discrete time-scale. Precisely, we give explicit 
formulas showing how behave the scale derivatives under change of scale. This is illustrated 
with the extension of Okamoto’s function, called scale and multiscale Okamoto’s function. 
We also discuss why there exists different continuous models associated with the same scale 
equation and we apply it on the Newton’s equation. The formalism developed in this paper 
goes beyond the time-scale calculus which is naturally included as a special case. 
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T 

'^elem 

Elem(T) 

T 

C(T,R) 
C(T,R) 

F 

A@F 

elem,A 

Oa 

Oa,N 
Oa,N,m 

Okamoto 

a{X, t) 
a{X) 


'^Tmo 

X-ext 

"V* 

■^ext 

X^ 

a 

P 

P 

V 

A 

Aoo 

V 

Voo 


(^) 


At 


□ 

at 


Notations 

Finite and discrete time scale 

Elementary time scale i.e. time-scale with two elements 

Decomposition of T in elementary time scale 

Scale sequence 

Space of discrete functions 

Space of scale functions 

A scale function 

Action of an operator A on a discrete function F 

Time scale obtained with the action of A 

Okamoto’s function 

Multiscale Okamoto’s function 

Multiscale Okamoto’s function of order m 

Set of multiscale Okamoto’s functions 

Scale range from T^q to with mo < mi 

Pointwise power-law scale regime in t 

Local power-law scale regime 

Global power-law scale regime 

Global scale regime 

Extension of the scale function X 

Reference scale function of Xg^t 

Asymptotic model of the scale function X 

Regular part of the asymptotic model of the scale function X 

Eorward jump 

Backward jump 

Eorward graininess function 

Backward graininess function 

Eorward derivative 

Asymptotic forward derivative 

Backward derivative 

Asymptotic forward derivative 

Gauchy A-integral 

Box derivative 
Asymptotic Box derivative 
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